Limits on the accuracy of isoelectronic gravity measurements at short separation due 

to patch potentials 
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In force sensing experiments intended to measure non-Newtonian gravitational signals electrostatic 
patch potentials can give rise to spurious forces, torques, and noise. Undesired patch-induced 
interactions can lead to systematic effects which limit accuracy, and noise can place lower limits 
on precision. In this paper we develop the theory for electrostatic patch effects on isoelectronic 
experiments, where their mean effect is nullified by design. We derive analytical expressions for 
the patch force and torque power spectrum to estimate the limitations introduced by patch-induced 
signals. 
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The presence of patch potentials on metallic surfaces 
has important implications in various experiments, in- 
cluding measurements of gravity on elementary particles 
PQ, heating in ion traps [5J |3], ionization of Rydberg 
atoms [3], precision tests of general relativity in space 
|B], measurements of the Casimir force [7HID]. and 
searches for hypothetical forces [TTlJ21| . Patches lead 
to interactions that limit the accuracy and precision of 
such experiments. To date, the theoretical modeling has 
focussed on quantifying the mean effect of patches on 
measured signals [22 26 , but an in-depth analysis of the 
patch force/torque power spectrum is lacking. Here, we 
develop such a theory and apply it to the analysis of two 
ongoing experiments that aim to measure new hypothet- 
ical forces in which the mean patch effect is cancelled out 
by construction, but where variations of the surface po- 
tential still detrimentally affect the ability to constrain 
new physics. 

These two experiments, one being performed at the 
University of Birmingham [27] and another at IUPUI 
|28j . use the isoelectronic technique. In this approach 
a flat material sample with inhomogeneous mass density 
is covered with a sufficiently thick metal coating to elec- 
trically screen the underlying materials. This has two im- 
portant consequences: Both (i) the Casimir interaction 
and (ii) the mean patch force between the coated sam- 
ple and a force probe are independent of the location of 
the latter at fixed separation distance. Force signals that 
depend on lateral displacements can arise not only from 
the sought after mass-dependent interactions between the 
probe and the underlying inhomogeneous mass density, 
but from spatially dependent patch force fluctuations as 
well. Hence, although the isoelectronic technique elimi- 
nates Casimir and mean patch forces, it is still possible 
that patch effects contribute to the measured signal. In 
order to quantify the magnitude of patch systematics, we 



compute the electrostatic torque and force power spec- 
trum for the above two experiments and show that, for 
reasonable estimations of parameters describing the sam- 
ple's surface potentials, patch effects can be an important 
limitation to accuracy. 

The Birmingham Experiment and Patch Torque: This 
experiment searches for violations of the inverse square 
law of gravity by measuring the torque acting between 
two parallel disks [27]. As shown schematically in Fig. 
la, two identical disks (radius R$ — 40 mm, thickness 
10 //m) are fabricated with a pattern of 2048 pairs of 
gold/copper spokes, and over-laid with a gold film (thick- 
ness 1 /Ltm, not shown in Fig. la) for isoelectronic measure- 
ments. The 2048 mass pairs are grouped into 16 sectors 
of alternating phase to provide a further modulation of 
the signal. In the schematic in Figure la, for clarity, 
only 36 mass pairs are shown in 4 sectors of alternat- 
ing phase. The lower disk, the test mass, is attached 
to a torsion pendulum with a superconducting suspen- 
sion and the upper disk, the source mass, is attached to 
a micro-positioner. The spacing between the gold sur- 
faces is nominally D = 13 /im. The experiment has been 
designed to achieve a sensitivity of 10~ 17 N m after an 
integration period of 1 day. The experimental procedure 
comprises the measurement of the torque as a function of 
the rotation angle, which should be angle periodic with 
a frequency determined by the mass density modulation. 

To compute the patch torque power spectrum we first 
obtain an expression for the patch torque as a function 
of angular displacement 4> between the disks. Given that 
D <C Rq, we approximate the disks by infinite parallel 
planes and neglect finite size effects. The electrostatic 
interaction energy as a function of <f> for given fixed patch 
potentials V" a (x) on the disks (a = 1,2) located at z a = 



0,D, is 
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and £q is the permittivity of 
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the xy— plane, k = |k 
free space. The matrix .R7 1 rotates vectors lying in the 
xy-plane (see supplementary information). The torque is 
t{4>) — d£(4>)/d(/). Notice that terms containing E do not 
contribute to the torque because they are independent of 
the relative angle. 

Without a precise knowledge of the surface potentials 
V a (x), we proceed by making some statistical arguments. 
For this purpose we adopt the "quasi- local" model of ran- 
dom surface potentials |25] . A given micro-realization of 
surface voltages in this model corresponds to a random 
tessellation of the surface with random and statistically 
independent voltage assignments to each patch domain. 
In the following we assume that the patch voltages fluctu- 
ate around zero, as a constant background potential does 
not contribute to the signal in the isoelectronic approach. 
Provided many elementary patch areas are contained on 
the plates we can replace the product of surface voltages 
V a (x)Vb(x') by the two-point voltage correlation func- 
tion, (V Q (x)V fc (x')) = <5 ab C aa (x,x') (ergodic hypothesis). 
The angled brackets denote ensemble averaging over all 
realizations of voltages and geometrical patch layouts. 
Consequently, (V 2 (x)Vi(x')) = and (t ((/>)) = 0. Be- 
cause the tesselation of the surface is determined ran- 
domly the voltage correlation exhibits translational and 
rotational invariance, i.e. C aa (x,x') = C aa (r), with r = 
|x — x'|. The explicit form of the quasi-local correlation 
function is given in Eq.(14) of [55]. For the special case 
when all patches have the same typical size £, and both 
disks have the same voltage correlation properties, this 
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is given by C ab (r) = ^^-[cos- 1 (|) - f (1- (f) 2 ) 1 / 2 ], 
where VJ- ms is the root-mean-square of the patch voltage. 
Unfortunately this correlator is inexpedient for the pur- 
poses of computation. Instead, we approximate it by a 
Gaussian 



C aa (x,x') « K 2 ms ex P (-4|x - xf/£ 2 ), 



(2) 



which captures the salient features of the quasi-local cor- 
relator and allows an almost analytical treatment (see 
inset of Fig. lb). 

The patch torque signal, although systematic, is ran- 
dom in the experiment. We quantify its variance by 
means of the patch torque power spectrum defined as 
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FIG. 1. (Color online) (a) Set-up of the Birmingham experi- 
ment to measure non-Newtonian gravity with a torsion pen- 
dulum using the isoelectronic technique. For clarity, the gold 
coating on top of the mass-density modulations of both disks 
is not shown. The underside of the top disk is also density 
modulated, (b) Root-mean-square amplitude of torque fluc- 
tuations as a function of typical patch size in the band rele- 
vant to gravity measurements in the experiment (m = 2048 
in Eq.Q). The solid line is for Vr ma = 45 mV (rms voltage 
for clean gold surfaces), and the dashed line is for V TnlB = 10 
mV (reduced rms voltage to mimic surface contamination). 
The gray horizontal line represents the design sensitivity of 
the experiment. The other parameters are D — 13 fira and 
Ro = 40 mm. The inset shows the quasi-local voltage corre- 
lation function (solid) and the Gaussian correlation function 
Q (dashed). 



Using the above correlation function we obtain 



Sr[m] = ^ 2 e 2 m 2 R 2 V^ 
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where J m (x) is the mth order Bessel function of the first 
kind [5S] (see supplemental information for details). The 
relevant value of the integer m is 2048 because this corre- 
sponds to the angular frequency of the mass-modulated 
signal. Using this expression we compute the expected 
torque fluctuations as a function of patch size and com- 
pare with the design sensitivity of the Birmingham ex- 
periment (see Fig. lb). For clean gold surfaces the rms 



patch torque fluctuations exceed the design sensitivity 
over nearly the entire range of physical patch sizes, imply- 
ing that patch effects dominate non-Newtonian signals. 
Surface contamination is known to lower the typical root- 
mean-square variations of surface patch potentials |30j . 
We mimic this effect with a reduced V rins in Fig. lb, 
which shows that the patch signal can be lowered by a 
factor of 20 when the rms potential variation is decreased 
by a factor of 4.5. This suggests that purposefully con- 
taminating the sample surfaces with adsorbates might 
beneficially suppress the patch torque power spectrum. 

The IUPUI Experiment and Patch Force: In this 
experiment a Au-coated sapphire sphere (radius R — 
150 /xm) is glued to a microtorsional oscillator (MTO) 
[2"8] , The sphere-MTO assembly is brought into close 
proximity to a rotating disk, as schematically shown in 
Fig. 2a. The apex of the sphere is at a distance L — 1 cm 
from the axis of rotation. The disk is the source mass 
formed by 1p (p = 200) sectors that subtend the same 
angle, which are alternatively made of Si and Au. All 
the sectors are covered by a uniform layer of Au (thick- 
ness 200 nm), thick enough to make the setup isoelec- 
tronic. The disk is set to rotate at an angular frequency 
fi = ui r /p, where uj r ca 2-7T x 300 Hz is the resonant fre- 
quency of the MTO-sphere assembly. The experiment 
is designed to achieve a sensitivity of 0.1 fN with an in- 
tegration time of 1000 s. The experimental procedure 
comprises measuring the normal force on the MTO at its 
resonant frequency ui r as a function of the disk-sphere 
separation D. 

The Fourier component of the patch force at the res- 
onance frequency of the MTO adds a systematic signal 
to the measurement. We estimate the magnitude of this 
effect by considering the patch force signal as a function 
of the rotation angle and describe the surface potentials 
using the approximate correlation function pi). For sim- 
plicity, we approximate the disk by an infinite plane, and 
treat the sphere's curvature using the proximity force ap- 
proximation (PFA). The latter assumes that the force be- 
tween non-planar objects is the sum of the forces between 
infinitesimal planar sections, and relates the sphere-plane 
force F to the plane-plane energy, £, as F = 2ttR £/A, 
where A is the area of the plate. In [21] the validity of 
the PFA for describing patch effects in the sphere-plane 
geometry was tested and verified for D <C R and for a 
variety of patch sizes (see Fig. 4 in this reference). The 
plane-plane energy is given by 



£((/>) = [ d 2 x [ d 2 x'V 1 (x.)A(x-x',D)V 2 (-x') 
Ja J a 



(5) 



where now x = R7 ■ (x + Lx) — Lx and x is a unit vector 
pointing in the ^-direction. The origin of the Cartesian 
coordinate system is centered on the point of closest sepa- 
ration between the sphere and the disk, so that the series 
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FIG. 2. (Color online) (a) Set-up of the IUPUI experiment 
to measure non-Newtonian forces with a torsional balance us- 
ing the isoelectronic technique. For clarity, the gold coating 
on top of the mass density modulated disk is not shown, (b) 
Expected force fluctuations in the band relevant to gravity sig- 
nals in the IUPUI experiment as a function of typical patch 
size for different separations D: 100 nm (solid), 150 nm (long- 
dashed), 200 nm (short-dashed) and 250 nm (dotted). K ms is 
set to 45 mV. The gray horizontal line demarcates the force 
sensitivity in the IUPUI experiment at 0.1 fN. The inset ap- 
proximates the effects of contamination which maintains the 
same convention as the main figure except Kms is set to 10 
mV. 



of transformations leading to x represents a rotation of 
the plane on an axis centered at —Lx. The dots indi- 
cate terms independent of the rotation angle that do not 
contribute to the relevant patch signal. Note that the 
mean value of the angle-dependent energy is zero since 

(y 1 (x)v^ 2 (x')) = o. 

The patch force power spectrum is defined as 



SpM 



dt' e lu,t (F(t + t')F(t)), 



(6) 



where T is the time period over which the measured sig- 
nal is averaged, and the time-dependence of the force 
arises from replacing </> with ilt in |5]). The force power 
spectrum at the resonance frequency of the MTO is then 
given by 



3 ff 2T/4 p- 

b F [u r \ - — s 



t/-fc" 



dk- 



A J sinh 2 kD 

x[j2(ky/A/^) - J p+1 (fc v ^r)J p _i(fe v /I7^r)]. (7) 



Note that it depends upon the plate area, which is ex- 
pected given that PFA can be viewed as an area average 
of the plane-plane force (see supplementary information) . 
The finite size of the sphere limits the sphere-plane in- 
teraction to an effective area A e g on the plane. Within 
the PFA we can estimate it by equating the sphere-plane 
force to the product of the plane-plane pressure P as 
F = PA cS . Using that E/A ~ DP, it follows that 
A e Q = 2-kRD. Therefore, the relevant area in ^f\ is 
A = 2n RD. An important observation is that the patch 
force variance depends upon the rotation frequency (re- 
call p = ui r /ft). Hence the influence of patches can be 
ruled out if one measures the variance of the force signal, 
assuming no other systematic signals are present, for two 
different rotation frequencies and their ratio is unity. For 
example, the ratio of Sp[u r ] for rotation frequencies uj r /p 
and co r /(p + 5) is 0.77 at D — 150 nm and for 100 nm 
patches. By taking the square root of Sf [w r ] we can es- 
timate the magnitude of the systematic patch force fluc- 
tuations at the frequency where the sought-after gravity 
signal should be present. In Fig. 2b the magnitude of the 
patch signal is computed for several sphere-plate separa- 
tions. We find that patch force signals are large enough to 
be measured for clean gold surfaces in the distance range 
150 — 250 nm and therefore may detrimentally limit con- 
straints on hypothetical new forces. The inset shows that 
contamination may prevent patches from influencing the 
IUPUI experiment altogether. 

A similar analysis can be performed for the previous 
isoelectronic experiment done at IUPUI in 2005 [15]. In 
this experiment the inhomogeneous source mass was de- 
posited on a MTO and the sphere was moved from re- 
gions of high to low mass density and back, following 
a periodic motion z(t) = D + SD cosoj z t normal to the 
MTO and x(t) = xq cosui x t parallel to its axis. The mod- 
ulation frequencies summed to the resonant frequency ui r 
of the MTO. The experimental noise was observed to be 
0.3 fN with an integration time of 1000 s. The experiment 
measured a non- vanishing signal (see Fig. 3 of [H]), that 
was not attributed to a new hypothetical force but to a 
differential Casimir force arising from a small variation 
in the height of the coated sample on the MTO. Even in 
the absence of such height variation, random (although 
systematic) patch forces could have contributed to the 
measured signal in Fig. 3 of [15]. For this set-up the 
patch force power spectrum at the MTO's resonance fre- 
quency is given by (see supplementary information) 
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where C[k] is the Fourier transform of pi) and, as before, 
A = 2irRD. In Fig. 3 we plot the magnitude of the 
expected patch force fluctuations as a function of patch 
size for different sphere-MTO distances. It is possible 
that the signal observed in [T^ is due to patches. The 
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FIG. 3. Expected force fluctuations in the band relevant to 
gravity signals for the 2005 IUPUI experiment [T5] as a func- 
tion of typical patch size for different separations D: from top 
to bottom 100 nm, 250 nm, and 500 nm. The gray horizontal 
line indicates the force sensitivity for the experiment. The 
rms patch potential is set to that of clean gold, V rms = 45 
mV. The inset shows the expected patch force fluctuations at 
150 nm (dashed line) and the black horizontal line represents 
the magnitude of the measured signal in Fig. 3 of [TS] at 
this particular distance (as in the main figure, the gray hori- 
zontal line is the force sensitivity). The constraints on forces 
reported in [15] imply that patches cannot have sizes exceed- 
ing 600 nm (corresponding to patch sizes right of the vertical 
line), otherwise they could result in patch forces larger than 
the total measured signal. Other parameters are: 8D = 10 
nm, xq = 75 /jm, and R = 50 fim. 



data of Fig. 3 together with Fig. 3 of [IS] can be used 
to constrain the possible patch sizes which could hypo- 
thetically produce the observed force. In order for our 
model to be compatible with the observations in [TS] the 
patches must be smaller than 600 nm for clean gold. 

In summary, we have developed the theory for patch ef- 
fects relevant to isoelectronic experiments aimed at mea- 
suring new hypothetical forces, in which the mean effects 
of patches are nullified by design, but whose fluctuations 
induce non-zero systematic patch force and torque sig- 
nals. We have found estimates for the root-mean-square 
amplitude of such fluctuations, and quantified limits on 
the accuracy of two proposed isoelectronic set-ups. The 
application of our analysis to a previous isoelectronic ex- 
periment suggests the possibility that part of the force 
signal could have been due to the presence of patches on 
the samples. Our analysis indicates that the systematic 
effects of patches in all these experiments decrease for 
contaminated samples. Although the patch model em- 
ployed in this work is inspired by physical arguments, the 
actual experimental patch distribution need not be pre- 
cisely described by this model. Consequently, our limits 
on accuracy to isoelectronic experiments are only quali- 
tative, a quantitative analysis requires the measurement 
of the patch distribution for the given experimental sam- 
ples. 
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